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This talk summarizes the quenched calculation of /s and /s a presented in [jlj. The heavy quark is simulated 
using an 0(1/M 2 ) improved NRQCD action, and the tadpole improved clover action is used for the light quarks. 
The axial current includes the 0(1/M) corrections in one-loop perturbation theory and the 0(1/M 2 ) ones at 
tree-level. We discuss the various systematic effects and consistency checks made in the calculation. 
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1. NRQCD ACTION AND OPERATOR 

The forward-backward symmetric evolution of 
the heavy quark propagator Gt using the NRQCD 
action |2] is given by 



Gt-i 
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where i?o = A^ 2 ' /2toq and <5 ± i? is defined in [Q. 
The gauge links are tadpole improved by divid- 
ing them by uq — 0.87779. The tree-level contri- 
bution exp(— iTigt) of the bare mass itiq to the 
heavy quark propagator is omitted from the cal- 
culation. The leading correction proportional to 
a s redefines ttiq: we take this into account per- 
turbatively when matching the lattice theory to 
the continuum. 

The current J, expanded in powers of 1/M, is 



J = 



m j(0) + J« + Jl_ J(2 

' 2M 2Af 



'/2 



+ ^ 8) +^ 1) + ^ 6) . (2) 



Vi9 



8M 



8M 



UP 



where we keep all operators of dimension three 
and four, and only the three dimension five oper- 
ators that appear at the tree level: 

J(0) = q lhli Q, 

= g(7-A)7574Q, 

*In collaboration with A. Ali Khan (OSU), S. Collins 
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q is the four component light spinor and Q is the 
nonrelativistic heavy quark spinor whose upper 
two components are zero. 

On transcribing J on to the lattice, the mixing 
induced by the hard cutoff has to be included. 
Dropping dimension five and higher operators, 
the lattice currents Jl can be written as 




' Coo a Coi a Co2 

a C10 Cn C12 

\a C20 C21 C22 
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In the final expression for J ||, every a -1 ap- 
pears multiplied by an M _1 ; all infrared diver- 
gences cancel; and the final logarithmic scale de- 
pendence is of the form logawg. Apart from 
these logarithms, all the coefficients are finite for 
1 < am,Q < 00. We shall highlight the contribu- 
tion of the different currents to Jb- 

2. LATTICE PARAMETERS 

The statistical sample consists of 102 16 3 x 48 
quenched lattices at (3 = 6/g 2 — 6.0. Heavy 
quarks are simulated with amPq = 1.6,2.0,2.7 with 
n = 2; and 4.0,7.0,10.0 with n = 1. Light quarks 
with k = 0.1369, 0.1375 and 0.13808 (brack- 
eting the strange quark) were simulated using 
the tadpole improved clover action. Details of 
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the simulations and fixing of the lattice param- 
eters are given in Using the data for m\ 
and uip versus k we get ki — 0.13917(9) and 
a -1 = 1.92(7) GeV. For the strange quark mass 
we use k s = 0.1376(1) from mx to quote the 
central value, and n s = 0.1372(3) from rriK* to 
estimate the error. Lastly, a reanalyses of the B 
spectrum gives aM° = 2.32(12) [Q, whereas in 
§ we had used aM b ° = 2.22(11). This revised 
aM®, however, leads to an insignificant change in 
estimates of decay constants. Estimates of fs e , 
which only require an interpolation in k and rriq , 
are more reliable than those for fs , which need a 
large extrapolation in k. 

3. METHOD OF ANALYSIS 

To estimate uncertainty in measuring the mass 
and amplitude of two point correlators, we ex- 
tract the decay constants in three ways: 

1. Individual local-smeared correlators (JiS) 
and the smeared-smeared correlators (SS) 
are fit using the form Aexp— Et. From 
these we get the contribution of each bare 
operator to the decay constant: fiVM = 
Aij \/ Ag- These ji\f~M are combined us- 
ing Eq. H to obtain the decay constant for 
each heavy-light combination. Finally, we 
extrapolate/interpolate in k and arriq. 

2. We combine the Ji according to Eq. || to ob- 
tain (JS). We fit this and the (SS) to ex- 
ponential forms Aexp— Et and obtain the 
decay constant fV~M at each combination 
of heavy and light quarks. These are then 
interpolated/extrapolated in k and arriq. 

3. We make simultaneous fits to each (JiS) 
and the (SS) over the same fit range 
to obtain the /jVM. We then ex- 
trapolate/interpolate these to the physical 
masses. Finally, we combine these using 
Eq. H to obtain Jb and fs, ■ 

We find that for am Q = 1.6, 2.0 and 2.7, all three 
methods of analysis give the same result, whereas 
for arriQ = 4.0 the agreement is marginal. Conse- 
quently, we use only the first three values of amPq 
to obtain our final results. 



Table 1 

Various contributions to fq ^/Mq versus am\ 



am Q 


j(0) 




j(3) 


J(4) 


j(5) 


1.6 
2.0 
2.7 
4.0 


.164(7) 
.167(9) 
.170(11) 
.174(11) 


-.025(1) 
-.021(1) 
-.016(1) 
-.011(1) 


-.0053(4) 
-.0036(2) 
-.0022(1) 
-.0011(1) 


.0036(2) 
.0022(1) 
.0011(1) 
.0005(0) 


-.0050(3) 
-.0033(2) 
-.0018(1) 
-.0008(1) 



Table 2 

Jq^/Mq extrapolated to ki and n s . Results are 
given for aq* ~ 1 and 7r. 



am° Q 




K 


s 


l/a 


n/a 


l/a 


it /a 


1.6 


0.114(6) 


0.121(6) 


0.136(4) 


0.144(4) 


2.0 


0.121(7) 


0.128(7) 


0.144(5) 


0.153(5) 


2.7 


0.127(9) 


0.135(9) 


0.153(6) 


0.164(6) 


4.0 


0.135(9) 


0.144(10) 


0.164(6) 


0.176(6) 



Estimates of fi^/M, multiplied by their tree- 
level coefficients and extrapolated to m, are pre- 
sented in Table 0. We expect that the bare 
matrix element of /2toq should be smaller 
than that of by a factor of 0(a/ arnq) ~ 
0(kqcD / to q)j an d qualitatively this is borne out 
by the data. Similarly, the bare matrix elements 
of J^/8{m° Q ) 2 , J^/8{m Q ) 2 and J^/A(m Q ) 2 
are smaller than that of by a factor of 

0(l/am%). 

4. RESULTS 

Table || presents /q^/Mq extrapolated to ki 
and k s for two choices of the lattice scale we be- 
lieve covers the range relevant to this calculation, 
aq* — 1, 7r. (The lattice coupling at these scales 
are av^/a) = 0.1557 and a v (l/a) = 0.2453.) 
For our final values we average the two estimates 
and interpolate to aAf 6 °: 

Ib = 147(ll)(i? a )(09)(6)MeV 
f Bm = 175(08)(i[ )(ll)(7)(15)MeV 

^ = 1.20(4)^). ( 5 ) 

JB 

The error estimates arc obtained as follows. To 
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determine the uncertainty due to setting the lat- 
tice scale, we repeat the calculation with a -1 = 
1.8 and 2 GeV 0. This gives a variation (1® 2 ) 
MeV in f B and (±J ) in f Bm . 

To estimate the perturbative errors we consider 
three issues, (i) The results at q* = 1 and 7r dif- 
fer by about 5%. (ii) The dimension five opera- 
tors have been included with tree-level coefficients 
and their mixing with lower dimensional oper- 
ators neglected. This, we estimate, introduces 
an error of 0{a s /(Ma) 2 ) ~ 4%. (hi) The ne- 
glected higher dimension operators are expected 
to be 0{K QCD /Mf < 1%. Altogether we assign 
a ~ 7% error coming from these sources. 

The remaining discretization errors are 
O(Aqcdo) 2 , which we estimate to be 4%. The 
last error in fs s and fg s / fs is due to the uncer- 
tainty in k s . We use the difference in k s obtained 
from K and K* to estimate this. 

An alternate approach is to neglect the con- 
tribution of the dimension five operators j( 3 >, 
j( 4 ) and j( 5 ) in the final result, and use them 
only to estimate errors. The reason is that 
the neglected mixing with J*- -* could introduce 
0(a s /(Ma) 2 ) errors, which can be larger than 
their 0{Kqcd/M) 2 contribution. Neglecting the 
dimension five operators (in which case the cal- 
culation is consistent to 0(1/M) and O(afM)) 
would give 

J B = 152(1 1)(± 8 12 ) (10) (6)MeV 

f Bs = 181(08)(i; )(12)(7)CtS)MeV. (6) 

Finally, we give a breakup of the various con- 
tributions. The tree-level estimate is Jb = 169 
MeV. Of this 195 MeV comes from -21 
MeV from jW, and -5 MeV from the l/M 2 cor- 
rections. The one-loop corrections reduce the es- 
timates by ~ 13%. On the other hand we find 
that the ratio /b,//b is insensitive to the various 
corrections. 

5. DISCUSSIONS 

In Table || we compare our results with other 
calculations. The formalisms used for incorporat- 
ing heavy quarks are NRQCD H, Fermilab Hf7|, 
Wilson/clover ||], and Wilson ||. Within errors 
the different approaches give consistent results for 



Table 3 

Comparison with other lattice determinations. 
The three groups of results use NRQCD, Fermi- 
lab and Wilson or Clover formalisms respectively. 
The errors are statistical and systematic (exclud- 
ing quenching) combined in quadrature. 





fs 


fB. 


/b, / Jb 


Ill 

|5| 


147(ll)(t£) 
162(7X1?*) 


175(8)(16) 
190(5)(tg) 


1.20(4)(^) 
1.18(3)(5) 


|S| 

1 


164(tiJ)(8) 
173(4) (12) 


185(^ 3 )(9) 
199(3) (14) 


1.13(H) 


121 


180(32) 
159(ll)(lf) 


205(15) (31) 
175(10)(1*}) 


1.14(8) 
l.ll(2)(±t) 



the quenched theory. Further improvement of the 
NRQCD estimates needs the full operator mixing 
matrix at 0(1/M 2 ), which in turn requires im- 
proving the light quark action to 0{a 2 ). A neces- 
sary check of our results is to show that they are 
independent of a. We are addressing these issues. 

A surprising feature of these data is that the 
ratio /b s //b) which is insensitive to most of the 
systematic errors, shows a difference between the 
NRQCD and the other formalisms. We do not un- 
derstand this. Possible issues to investigate are (i) 
the chiral extrapolation and (ii) the dependence 
on the heavy quark mass as this is handled dif- 
ferently in the different approaches. 
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